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Abstract 

We address the impact of sterile fermions on the lepton flavor violating decays of 
quarkonia as well as of the Z boson. We compute the relevant Wilson coefficients 
and show that the 'B{V —> where V — y;A) , Z can be significantly 

enhanced in the case of large sterile fermion masses and a non-negligible active- 
sterile mixing. We illustrate that feature in a specific minimal realization of the 
inverse seesaw mechanism, known as (2,3)-ISS, and in an effective model in which 
the presence of nonstandard sterile fermions is parameterized by means of one heavy 
sterile (Majorana) neutrino. 


PACS: 14.60.St, 14.60.Pq, 13.20.Gd, 13.38.Dg. 



1 Introduction 


So far no signal of new physics has been observed but its search is important in order 
to understand how to enlarge the Standard Model (SM) to solve both the hierarchy and 
the flavor problems. One of the most signihcant observations requiring us to go beyond 
the Standard Model is the assessment that neutrinos are massive and that they mix (^. 
Possible SM extensions aiming at incorporating massive neutrinos give rise to interesting 
collider signatures and open the door to new phenomena such as lepton flavor violating 
(LFV) decays. 

Currently, the search for manifestations of LFV constitutes a goal of several experi¬ 
mental facilities dedicated to rare lepton decays, such as ia —t and i and to 

the neutrinoless /i — e conversion in muonic atoms. One of the most stringent bounds from 
these searches is the one derived by the MEG Collaboration, B(/i —?• ey) < 5.7 x 10“^^ [^, 
which is expected to be improved to a planned sensitivity of 6 x 10“^^ [^. Moreover, the 
bound B(/r ^ eee) < 1.0 x 10“^^, set by the SINDRUM experiment [^, is expected to be 
improved by the Mu3e experiment where a sensitivity ~ 10“^® is planned [^. Limits on 
the T radiative decays and the three-body decays of r appear to be less stringent 
right now, but are likely to be improved at Belle II j^, where the search for LFV decays of 
the B-meson will be made too j^. The most promising developments regarding LFV are 
those related to the /r — e conversion in nuclei. The present bound for the /U“Ti ^ e“Ti 
conversion rate is 4.3 x 10“^^ (l0|, and the planned sensitivity is ~ 10“^® 11 . Similar is 


the case for gold and aluminum (12,13 


Searches for LFV are also conducted in high-energy experiments and a first bound on 
the Higgs boson LFV decay h —?• /rr has been reported by the CMS Collaboration 
The LHCb Collaboration, instead, reported the bound B(r 3 fi) < 8.0 x 10 


is likely to be improved in the near future 16 . Notice also that they already improved the 


14 


15 , which 


bounds on B(5(5) —?• e/r) by an order of magnitude (17|. 


In this work we will focus on the indirect probes of new physics through the LFV 
processes of neutral vector bosons, namely V — >■ ia ^/ 3 , with e {e, /r,T}, and V G 
{0,Z}, where stands for J/'ip and its radial excitations, and similarly for 
Most of the research in this direction reported so far is related to the Z —>■ 
decay modes. More speciflcally, the experimental bounds, obtained at LEP are found to 
be B{Z < 1.7 x 10“® Q, B(Z ^ < 1.2 x 10“^ @[1^, and B{Z 

e^T^) < 9.8 X 10 ® (l^[^ . One of these bounds has been improved at LHC, namely 

B(Z —>■ < 7.5 X 10“^ [^. On the theory side, the Z decays have been analyzed 

in the extensions of the SM involving additional massive and sterile neutrinos that could 

mix with the standard (active) ones and thus give rise to the LEV decay rates 22-24 . A 


similar approach has been also adopted in Ref. 25 , in the perspective of a Tera-Z factory 
ECC-ee 26 for which a targeted sensitivity is expected to be B(Z ^ 


10 


-13 


27 


Lepton flavor conserving decays of quarkonia have been measured to a high accuracy 
which can actually be used to fix the hadronic parameters (decay constants). Otherwise, 
one can use the results of numerical simulations of QCD on the lattice, which are nowadays 
accurate as well 28-31 . The experimentally established bounds for the simplest LEV 
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decays of quarkonia are 


16 


B(0 ^ e^) < 2.0 X 10“®, 



B(J/^ ^ eii) < 1.6 X 10"^ 
B(J/^ ^ ^r) < 2.0 X 10“®, 

B(T ^ ht ) < 6.0 X 10■^ 


B(J/^ ^ er) < 8.3 x 10“®, 

gm 


35 


B(T(25) ^ er) < 8.3 x 10"®, 
B(T(3,S) ^ er) < 4.2 x 10“®, 


B(T(25) ^ iir) < 2.0 x 10“®, 
B(T(3,S) ^ fir) < 3.1 x 10“®, 


( 3 ^ 


where each mode is to be understood as B(l^ —>■ = B(l/ —>■ + B(l/ —>■ 

Despite the appreciable experimental work on the latter observables, only a few the¬ 
oretical studies have been carried out so far. The authors of Ref. 37 applied a vector 


meson dominance approximation to /r —?• 3e and expressed the width of the latter process, 
r(/r —>■ eee) = r(/r —>• Ve)T{V —>■ ee). Since the values of r(R ee) are very well known 
experimentally (^, the experimental bound on r(/r ^ 3e) is then used to obtain an upper 
bound on the phenomenological coupling gv^e) which is then converted to an upper bound 
on r(l/ fie). A similar approach has been used in Ref. where instead of /i —?• eee, 
the authors considered the fi — e conversion in nuclei (N), which they described in terms 
of a product of couplings gvue and gvNN- The latter could be extracted from the exper¬ 
imentally measured r(l/ —>■ pp), and with that knowledge the experimental upper bound 
on R(/iTi —7- eTi) results in an upper bound on r(R —?• /re). A more dynamical approach 
in modeling the V processes has been made in a supersymmetric extension of the 

SM with type I seesaw |M| . 

Sterile fermions were proposed in various neutrino mass generation mechanisms, but the 


interest in their properties was further motivated by the reactor/accelerator anomalies 40 - 


43 


a possibility to offer a warm dark matter candidate |44| - [46| , and by indications from 
the large scale structure formation (47| - |4^ . 

Incorporating neutrino oscillations (masses and mixing (^) into the SM implies that 
the charged current is modified to 


-£ec = + C.C. , 


( 1 ) 


{/ being the leptonic mixing matrix, a the flavor of a charged lepton, and i = 1,... ,71^, 
denotes a physical neutrino state. If one assumes that only three massive neutrinos 
are present, the matrix U corresponds to the unitary Pontecorvo-Maki-Nakagawa-Sakata 
(PMNS) matrix. In that situation the GIM mechanism makes the decay rates B(l/ —>■ 
completely negligible, < 10“^°. That feature, however, can be drastically changed in the 
presence of a non-negligible mixing with heavy sterile fermions. In what follows we will 
consider such situations, derive analytical expressions for B(R ^a^/s), and discuss a spe¬ 
cific realization of the inverse seesaw mechanism, known as (2,3)-ISS [^. We will also 
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discuss a simplified model in which the effect of the heavy sterile neutrinos is described 
by one effective sterile neutrino state with non-negligible mixing with active neutrinos. 
Despite several differences, our approach is similar to the one discussed in Ref. where 
the SM has been extended by new, heavy, Dirac neutrinos, singlets under SU (2) x C(l). 
and applied to a number of low energy decay processes. Our sterile neutrinos are Majorana 
and we apply the approach to the leptonic decays of quarkonia for the first time. 

The remainder of this paper is organized as follows: In Sec. we formulate the problem 
in terms of a low energy effective theory of a larger theory which contains heavy sterile 
neutrinos, we derive expression for B(R and compute the Wilson coefficients. In 

Sec. 1^ we briefly describe the specific models with sterile neutrinos which are used in this 
paper to produce our results presented in Sec. We finally conclude in Sec. 


2 LFV decay of Quarkonia - Effective Theory 

In this section we formulate a low energy effective theory of the LFV decays of quarkonia 
of type V and express the decay amplitude in terms of the quarkonium decay 

constants and the corresponding Wilson coefficients. The latter are then computed in 
the extensions of the SM which include the heavy sterile neutrinos. We also derive the 
expression relevant to r(Z 


2.1 Effective Hamiltonian 


Keeping in mind the fact that we are extending the SM by adding sterile fermions, with¬ 
out touching the gauge sector of the theory, the decays of vector quarkonia, V{q) —?• 
^aip) “ P); only occur through the photon and the Z-boson exchange at tree 
level. In the lepton flavor conserving processes the Z-exchange terms are very small with 
respect to those arising from the electromagnetic interaction and are usually neglected. 
The generic effective Hamiltonian can be written as 


6^02 _ 

'^eS = QIiJ'Q ■ ^0 


Jf 


Cvl'J^Pl + Cvr^^Pr H- (CrPr + CrPl) 


Tn\Y 


h, ( 2 ) 


where Qq is the electric charge of the quark Q, my is the mass of quarkonium V which is 
dominated by the valence quark configuration QQRy r^l,r are the Wilson coefficients, 
p is the momentum of one of the outgoing leptons, and Pr/r = |(1 TTs)- Contributions to 
the scalar (left and right) terms are suppressed by ma,/ 3 /mw, where ma ,/3 are the charged 
lepton masses. In this section we will keep such terms so that our expressions can be useful 
to approaches in which the scalar bosons are taken in consideration. For our phenomeno¬ 
logical discussion, however, it is worth emphasizing that CRyyR <C Cvr. 

^In this work, due to the tension between the most recent Planck results on extra light neutrinos 
(relics) and the reactor/accelerator anomalies, we will consider the effect of (heavier) sterile neutrinos not 
contributing as light relativistic degrees of freedom . We will require our models to be compatible with 
current experimental data and constraints and to fulfill the so-called perturbative unitary condition which 
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Figure 1: Diagrams contributing the LFV decay of quarkonia V The blob in the first diagram 

is related to the penguin loop that generates the LFV, and the box diagram is particularly important to 
be included in the case of ^ ^a^/3 because of Vth ^ 1 and of the top quark mass, making the box 
diagram contribution to the Wilson coefficient significant. 


Without entering the details of calculation it is easy to verify that the only relevant 
diagrams are those shown in Fig. and therefore the structure of the Wilson coefficients 
Ci reads, 


|2 ^ ^ (o\ 

sm^Owm^QQ^ ^ ’ 

where (7/’^ are the contributions arising from the photon and the Z-boson exchange, while 
comes from the box diagram that involves the Cabibbo-Kobayashi-Maskawa coupling 
0 In the above expressions gy — — Qgsm^Ow The blob in the diagram shown 

in Fig. [^stands for the lepton loop diagrams that may contain one or two nentrino states 
and which, in the extensions of the SM involving a heavy neutrino sector, will give rise 
to the LFV decay due to the effect of mixing which is parametrized by the matrix U [see 
Eq. 01 - Separate contributions coming from different diagrams can be further reduced by 
factoring out the neutrino mixing matrix elements, namely 


a = c? + cf 


1 


m 


V 


9v 


sin"^ 9w cos^ 9w ~ "^1 Cq 




= ^u, 


riu 


13k 


KkC] 


,Box;A; 


and Cf = 


k=l 


k=l 


rijy n 

k=l j=l 


(3k^ 


(4) 


where we see that the term involving two neutrino eigenstates appears only in the Z 
coefficient because it is related to the vertex It is worth emphasizing that the 

tensor structure in Eq. (j^ can be easily obtained from the coefficients Cl,r by applying the 


52 


puts a strong constraint on the models for the very heavy sterile fermion(s) 

^We remind the reader that the ground vector meson cc, bb states are"^ T/'b, T, respectively, and 
the corresponding charges are Qs^b = — 1/3 and Qc = 2/3. 

^The box diagram contribution to F ^ in the case of F = T is dominated by the top quark 
1); for F = 'ijj it is negligible because the contribution of the b quark is Cabibbo suppressed 
{\Vch\ — 0.004) while the Cabibbo allowed one {\Vcs\ — 0.99) is suppressed by the strange quark mass; 
for F = 0, the contributions of the charm and top quarks are comparable but overall smaller than in the 
T ^ case. 
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Gordon identity. Such contributions are 1/mw suppressed, and thus completely negligible, 
which is why we do not give explicit expressions for these coefficients. 

Using the effective Hamiltonian Q and parameterizing the hadronic matrix as 

{^\QlixQ\V{q, cr)) = fvmvel , (5) 


where fy is the decay constant of a quarkonium V with momentum q and in a polarization 
state cr, we can write the decay rate as. 


r(V' ^ Cit) = 


Srriy 


fv 


with 

and 
4>c = 


\{a\ b\ c") = [a" - (6 - cY][a^ - (6 + c)% 


my ' 


2ir..2 


( 6 ) 


( 7 ) 




m 


V 


tr 


pH 


pH 


— f + nip) ■ {CyL^^PL + CyR^^Pn + Cl - Pl + Cr - Pr) 


TTiw 

y 


TTlyy 


m«) • {C*yLYPL + C*y^YPR + CI^Pr + C*^^Pl) 


TTlyy 


TTlyy 


( 8 ) 


which gives 


<Pc = 


(ruy -m\- + \CrY) - 4Re(G2Gij)m„m;3 

+ 4mvi/Re(C'2(C'yLm^ + CyRnia) + C*^{CyL'ma + CyRTtlj^)) 


2\2 


+ ^m^{\CyL\ + \CyR\Y 2my - my{m^^ + mj) - (m„ - mj) 


(9) 


+ 48m^mymam;3Re(G2^iGvii)|. 


As we mentioned above, we consider in our framework CyR ^ CyR^R^i, and therefore we 
can write 


r(U^£„£^)— fyCpmyr\CyL\ ' {my,m^,mYjY 2m? 2m^ 

( 10 ) 


V 


where \{a?,l?,C‘) is given in Eq. Q. In this last expression we also used r(U = 

Besides quarkonia we will also revisit the issue of adding extra species of sterile neutrinos 
to the decay of Z In that case the effective Hamiltonian can be written as 


^eff — 


9 


2 cos 9^? 


DyLl^Ph + DyRl^PR + DrPr + DrPr 




( 11 ) 
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Figure 2: Vertex diagrams contributing the LFV decays. 


where the Wilson coefficients are now denoted by D* and take the form 


d, = J 2 + E E . 

k=l j=l 


( 12 ) 


k=l 


The decay rate in the similar limit, Dvl ^ Dvr,r,l, reads 


-p//7 «-«+\ 8y/2 |2\l/2/ 2 2 

r(^ = iV—— on \^vl\ a ' {mz,m^,m 

^ oTTUlz COS^ U\^ 


^ + m}) _ - 1^1) 


2 \ 2 -\ 


2m| 


2m% 


(13) 


2.2 Wilson coefficients 


Concerning the computation of the Wilson coefficients we stress again that our results are 
obtained in a theory in which the Standard Model is extended to include extra species 
of sterile fermions, without changing the gauge sector. The origin of the leptonic mixing 
matrix U is model dependent and in order to be able to do a phenomenological analysis, 
we will have to adopt a specific model which will be discussed in the next section. 

The blob in the diagram shown in Fig. [T] stands for a series of diagrams such as those 
displayed in Fig. All of them, including the box diagram in Fig. have been computed 
in the Feynman gauge and the results are collected in Appendix A. Here we focus on the 
most important contributions in the case of large masses of sterile (Majorana) neutrinos. 
Contributions to the Wilson coefficients coming from vertex diagrams can be divided into 
two pieces: those involving only one neutrino in the loop, Cy2{xi), where Xi = 
and those with two neutrinos in the loop, Cy]^{xi,Xj). In the limit of large values of 
Xij ^ 1, we find the following behavior 




Cvd^i) -^ + finite term + 0{l/xi) logXj, 




327r2 

i»i a 
-> 


^ Cii Xi log X 


{2xi + 3 - dloga:,) + a:* ( loga:* - ^ ) |> + 


(14) 


To illustrate the relative contribution of the different diagrams we fix the values of the 
coefficients Cij = ^ U^^Uaj = 10“^, and plot \CvL{xi) — Cvl{^)\ and \Cyj^{xi,2 Xi) — 

OL = e^lJL^T 
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C'y^(0,0)| for the case of T —>■ /ur, cf. Fig. 10 We see that only for very large masses 
the diagrams with two neutrinos in the loop become more important than those with one 
neutrino state. We should stress that each contribution to CvLi^i), i.e. and 

Cyhixi), scales as logXj for large values of Xi, except for Cyj^{xi) which goes to a constant 
in the same limit. That can also be seen in Fig. where in the left panel we show the 
dependence of the total C'vi(xj) on Xi and in the right panel we show Cyj^{xi) and its 
dependence on the mass of the initial decaying meson, 0, and T. The contribution 

of sterile neutrinos to the LFV decay of T is larger than the one to lighter mesons, since 
the Wilson coefficients are also proportional to the mass of the initial particle. 




Figure 3: In the left panel are shown CvL{xi) and Cyi{xi, xj), for xj = 2xi, as functions of rrii = mw\/^, 
the mass of the heavy sterile neutrino propagating in the loops. For illustration purpose, the couplings 
Cij were fixed to a common value, 10“®, and the example corresponds to the T —> /rr decay. Right panel: 
Cybixi) is plotted as a function of ruj for the case of T —> e/r in three specific cases V G {<p, J/ip, T}. In 
both cases the value of functions at Xij = 0 have been subtracted away. 


Before closing this section we should reiterate that our Wilson coefficients have been 
computed in the Feynman gauge. Since all divergencies cancel out, our results are finite 
and gauge invariant, as was already observed in Refs. |22|-[^. 


3 SM in the presence of sterile fermions 

With the expressions derived above, we now have to specify a model for lepton mixing 
(couplings) Uci in the presence of heavy sterile neutrinos propagating in the loops. We opt 
for a minimal realization of the inverse seesaw mechanism for the generation of neutrino 
masses, which is nowadays rather well constrained by the available experimental data. 

^Due to the unitarity of the mixing matrix t/, the terms in the Wilson coefficients that do not depend 
on neutrino masses give a vanishing contribution after summing over all neutrino states. We thus subtract 
the constant terms in the plots in order to better appreciate the dependence on the neutrino masses. Notice 
also that Cij = 10“^ is in agreement with all constraints discussed in the text when the neutrino masses 
are below 0(100) TeV. 
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Furthermore, we will use a parametric model containing one effective sterile neutrino, 
which essentially mimics the behavior at low energy scales of mechanisms involving heavy 
sterile fermions. 


3.1 The (2,3)-inverse seesaw realization 

Among many possible realizations of accounting for massive neutrinos, the inverse seesaw 
mechanism (fSS) offers the possibility of accommodating the smallness of the active 
neutrino masses for a comparatively low seesaw scale, but still with natural 0{1) Yukawa 
couplings, which renders this scenario phenomenologically appealing. Indeed, depending 
on their masses and mixing with active neutrinos, the new states can be produced in 
collider and/or low energy experiments, and their contribution to physical processes can 
be sizable. ISS, embedded in the SM, results in a mass term for neutrinos of the form 

->Cmass = ^nlCMm + h.C. , (15) 

where C = is the charge conjugation matrix and Here i'L,a, « = 

e, /r, r denotes the active (left-handed) neutrino states of the SM, while ^R,i (* ~ h 
and Sj (j = 1, are right-handed neutrino fields and additional fermionic gauge singlets, 
respectively. The neutrino mass matrix M then has the form 


/ 0 d 0\ 

M= \ (f 0 n , (16) 

y 0 nF J 


where d, n, [i are complex matrices. 

The Dirac mass matrix d arises from the Yukawa couplings to the SM Higgs boson, H = 

ia^H, 


Y 





(17) 


while the matrix ii^ instead, contains the Majorana mass terms for the sterile fermions 
Sj. By assigning a leptonic charge L = +1 to both and s, one makes sure that the 
off diagonal terms are lepton number conserving, while s'^Cs violates the lepton number 
by two units. Furthermore, the interesting feature of this seesaw realization is that the 
entries of /i can be made small in order to accommodate for the 0(eV) masses of active 
neutrinos, with large Yukawa couplings. This is not in conflict with naturalness since the 
lepton number is restored in the limit of ^ 0. 

Concerning the additional sterile states ur and s, since up to now there is no direct 
evidence for their existence and because they do not contribute to anomalies, their number 


^It is in general possible to consider also a nonzero value for the central entry of the matrix (16), with 
elements at a mass scale similar to the one of r. These parameters, however, only affect neutrino masses 
and mixing at loop level 


55 


, which is why we do not consider them here. 

^In this work we consider configurations in which the entries in the above matrices fulfill a naturalness 
criterion, |/i| ^ \d\ < \n\ 
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is unknown. In Ref. 50 it was shown that it is possible to construct several minimal distinct 


realizations of ISS, each reproducing the correct neutrino mass spectrum and satisfying 
all phenomenological constraints. More specifically, it was shown that, depending on the 
number of additional fields, the neutrino mass spectrum obtained for each ISS realization is 
characterized by either two or three mass scales, one corresponding to rrii, ~ fi (P/n^ (light 
neutrino masses), one corresponding to the heavy mass eigenstates [the mass scale of the 
matrix n of Eq. (16)], and finally an intermediate scale ~ fi, only present if This 


allows us to identify two truly minimal ISS realizations that comply with all experimental 
bounds, namely the (2,2)-ISS model, which corresponds to the SM extended by two right- 
handed (RH) neutrinos and two additional sterile states, leading to a three-flavor mixing 
scheme, and the (2,3)-ISS realization, where the SM is extended by two RH neutrinos 
and three sterile states leading to a 3-|-l-niixing scheme. Interestingly, the lightest sterile 
neutrino with a mass around eV in the (2,3)-ISS can be used to explain the short baseline 
(reactor/accelerator) anomaly |40|-|4^ if its mass lies around eV, or to provide a dark matter 


candidate if the lightest sterile state were in the keV range 56 


3.2 A model with one effective sterile fermion 

Since the generic idea of obtaining a significant contribution to our observables applies to 
any model in which the active neutrinos have sizable mixing with some additional singlet 
states (sterile fermions), we can use an effective model with three light active neutrinos 
plus one extra sterile neutrino. 

The introduction of this extra state implies three new active-sterile mixing angles 
(^ 14 ,^ 24 ,^ 34 ), two extra Dirac CP violating phases (^ 14 ,^ 34 ) and one additional Majo- 
rana phase ( 04 i). The lepton mixing matrix is then a product of six rotations times the 
Major ana phases, namely 

u = ^34(^34; <^34) • -^24(^24) • -Ri4(^14; ^14) • -R23 • -Ri3 • Ri 2 • diag(021, 031, 04l) 

= ^34(^34, (^34) • -^24(^24) • -^14(^14, (^14) • DpMNS • diag(021, 031 , 04 l) , ( 18 ) 

where the rotation matrices R 34 , R 24 , -R 14 can be defined as: 


-R34 — 


Roa — 


Ria — 


/ 

1 

0 

0 



0 

\ 


0 

1 

0 



0 



0 

0 cosd 

34 

sind 34 • e 


V 

0 

0 — sin6* 

34 

. 01(^34 


COsd34 

/ 

/ 

1 

0 

0 

0 

\ 




0 

COsd24 

0 

sind; 

24 




0 

0 

1 

0 




V 

0 

— sind24 

0 

COS0; 

24 / 





COSdi4 


0 0 

sin 6 *i 4 • e 

\ 



0 


1 0 


0 




0 


0 1 


0 



— sin0i4 • 

0 0 


COsdi 4 

J 


(19) 
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In the framework of the SM extended by sterile fermion states, which have a nonvanish¬ 
ing mixing with active neutrinos, the Lagrangian describing the leptonic charged currents 
becomes 

-£ec = + c.c., (20) 

where i = 1,... ,nt, denotes the physical neutrino states, and a = e,fi,T are the flavors 
of the charged leptons. In the case of the SM with three neutrino generations, U is the 
PMNS matrix, while in the case of > 4, the 3x3 submatrix (t/pMNs) is not unitary 
anymore and one can parameterize it as 


Upum —t b^PMNS — (1 “ i?) b^PMNS ; 


( 21 ) 


where 77 is a matrix that accounts for the deviation of I7pmns from unitarity 


57 


58 


due 


to the presence of extra fermion states. Many observables are sensitive to the active-sterile 
mixing and their current experimental values can be used to constrain the 77 matrix 


59 


In order to express the deviation from unitarity in terms of a single parameter, we 
define 


77 — 1 — I det t/pMNsI) ( 22 ) 

which, in the case of the extension of the SM by only one sterile fermion and in terms of 
the mixing angles defined above, reads 

77 = 1 — I cos 614 cos 624 cos 634 1. (23) 

4 Results and discussion 

In this section we present and discuss our results. 

Since the Wilson coefficients of the processes discussed here are proportional to the 
mass of the decaying particle, it is quite obvious that the most significant enhancement 
of B(f/ ia^/3) will occur for f/ = T and its radial excitations. For this reason we will 
present plots of our results for this decay channel. Plots for other channels are completely 
similar which is why we do not display them. Before we discuss the impact of the active- 
sterile neutrino mixing on the LFV decay rates further, we first specify the constraints on 
parameters of both of our models. 

In Fig. 1^ (left panel), we plot the dependence of 77 with respect to the mass of the effective 
sterile neutrino 777 . 4 . Gray points in that plot are obtained by varying the mass of the lightest 
neutrino, G (10“^^, 1) eV, and by imposing the following constraints: (i) Neutrino data 
(masses and mixing angles) respect the normal hierarchy, with = 7 . 5 ( 2 ) x 10“^ eV, 

and A?? 7 .|]^ = 2.46(5) x 10“^ eV ID We checked to see that our final results do not 
change in any significant manner if the inverse hierarchy is adopted. Furthermore, we 
vary the three mixing angles with the fourth neutrino by assuming 9 i 4 G (0,27r], while 
keeping the other three mixing angles to their best-fit values, namely sin^ ^12 = 0.30(1), 
sin2d23 = 0.47(4), sin^^is = 0.022(1) ID (ii) The selected points satisfy the upper bound 
B(/i —7- ey) < 5.7 x 10 “^^ [^. (in) The results for = r(7r ehe)/r(7r ^ Rk, 
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0.01 
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Figure 4: 77 parameter, which parametrizes the size of mixing between the active and heavy sterile states, 
is plotted vs the mass of the heavy sterile state. The gray points (left panel) correspond to solutions 
complying with all experimental data and constraints discussed in the text except for perturbative unitary 
condition ([24|), which we then applied to obtain the region of dark-blue points. In the case of the (2,3)-ISS 


model (right-panel), we further imposed constraints of Ref. 59 on the matrix 77 , as well as the bound 
B{ii ^ eee) < 10 “^^, resulting in the bright-blue region of points. 


r(VF ^ and T{Z invisible), remain consistent with experimental findings. We see 
that for all (heavy) sterile neutrino masses the unitarity breaking parameter is 77 < 0.005. 
That parameter space is not compatible with the perturbative unitarity requirement, which 
for 7714 ^ '^w translates into 


GFTn\ 

\/27r 




< 1 . 


(24) 


The resulting region, i.e. the one that satisfies constraints (i), (ii), (iii) and Eq. (24), is 


depicted by bine points (the dark region) in Fig. where we see that the parameter r] is 
indeed diminishing with the increase of the heavy sterile mass 771 . 4 . In other words, the de¬ 
coupling of a very heavy sterile neutrino entails the unitarity of the 3x3 submatrix Vpmns- 
Decoupling from active neutrinos for very large masses was also explicitly emphasized in 
Ref. [^. We shonld mention that, besides the above constraints, we also implemented 
the constraint coming from B(/i ^ eee) < 10“^^ [^, but it turns ont that the present 
experimental bound does not bring any additional improvement. 


By imposing the constraints (i) and Eq. (24) on the (2,3)-ISS model, we get a similar 
region of allowed (blue) points in the right panel of Fig. A notable difference with respect 
to the sitnation with one effective sterile neutrino is that the region of very small mixing 
angles is excluded due to relations between the active neutrino masses and the active-sterile 
neutrino mixing, cf. Ref. [^. For very heavy 7775 , on the other hand, the range of allowed 
7 ]’s shrinks and eventually vanishes with 7775 00 . Furthermore, we use the results of 


^To write it in the form given in Eq. (24), we replaced aw = 5 ^/( 47 r) = \/2GFrti^ jn. 

®We recall that, in the (2,3)-ISS modeI~m 4 stands for the mass of the light sterile state whose impact 
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Ref. 59 which are derived in the minimal unitarity violation scheme in which the heavy 


sterile neutrino fields are integrated out, and therefore the observables computed in that 
scheme are functions of the deviation of PMNS matrix from unitarity only [^. We adapt 
and apply them to our (2,3)-ISS model and get a region of the bright-blue points, as shown 
in the right panel of Fig. To further constrain the parameter space we find it useful 
to account for the experimental bound on B(/r eee) < 1 x 10“^^, as is discussed in 
Refs. 23,60,62 . This latter constraint appears to be superfluous in most of the parameter 


space, once the constraints of Eq. (24) and Ref. |59) are taken into account, except in the 


range 10 TeV < < 100 TeV, where the bound B(/i eee) < 1 x 10“^^ restricts the 

parameter space relevant to B(l/ —>■ e/r). 

We also mention that we attempted implementing the constraints coming from various 
laboratory experiments, summarized in Ref. [^, but since those results only impact the 
region of relatively small sterile neutrino masses [m^ < 100 GeV), they are of no relevance 
to the present study. 


After having completed the discussion on several constraints, we present our results 
for branching fractions B(l/ ^ /ur) depending on the mass of heavy sterile neutrino(s). 
In Fig. [^we plot our results for R = T and V = Z, for which the enhancement is more 
pronounced. Other cases of V result in similar shapes but the upper bound becomes lower. 
In Table we collect our results for three values of the heavy sterile neutrino(s) mass. 

To better appreciate the enhancement of the LFV decay rates shown in Fig. we 
emphasize that both of them are B(R —>■ /it) < 10“™ in the absence of heavy sterile neu¬ 
trinos. Current experimental bounds in both cases are shown by dashed lines. Since those 
bounds are expected to improve in the near future, a possibility of seeing the LFV modes 
discussed in this paper might become realistic. Conversely, an observation of the LFV 
modes V —>■ ia^/ 3 , with branching fractions significantly larger than the bounds presented 
in Table would be a way to disfavor many of the models containing heavy sterile neutrinos 
as being the unique source of lepton flavor violation. In obtaining the bounds presented 
in Table [T] we used masses and decay constants fisted in Appendix B. In presenting our 
results (the upper bounds) for lepton flavor violating modes, we used the parameters from 
Ref. which were determined at 90% C.L. For that reason, we treated all other input 
data to 2 cr as well. Therefore, our final results in Table are also obtained at 2 cr level. 

Finally, we compare in Table our upper bounds for the modes for which we could find 
predictions in the literature. 

5 Conclusions 

In this paper we discussed the enhancement of the LFV decays of flavorless vector bosons, 
V with V G {0, %}, induced by a mixing between the active and sterile 

neutrinos. The enhancement grows with the mass of the heavy sterile neutrino(s), as can 

on the decays discussed here is negligible [as seen from Eq. (|l4|)], while ms > can be large and is 
important for B(T —)■ CC)- 
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Figure 5: B(T ^ /ir) and B(Z ^ jir) are shown as functions of the heavy sterile neutrino(s) mass, and in 
both models considered in this paper. The points are selected in such a way that the models are consistent 
with the constraints discussed in the text and shown in Fig. Dashed horizontal lines correspond to the 
current experimental upper bounds for these decay rates. Notice again that the mass of the heavy sterile 
neutrino is denoted as 7714 in the effective model, and 7775 in the (2,3)-ISS model because the latter contains 
a lighter sterile neutrino state, the impact of which is negligible on the decay modes discussed here. 
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Table 1: Upper bound on B(U —> laf-p) for three values of the mass 7714 ^ 5 . The numbers in the three 
columns referring to 7774 are obtained by using the effective model discussed in the text, while the other 
three, referring to 7775 , are results of the (2,3)-ISS model (also discussed in the text). 


be seen from the mass dependence of the Wilson coefficients that we explicitly calculated. 
We find that the most significant diagram that gives rise to the LFV decay amplitudes is 
the one coming from the Zuu vertex, which suggests a steady growth of the decay rate 
with the mass of the sterile neutrino(s). In the physical amplitude, however, the region of 
very large mass of the sterile neutrino(s) is suppressed as the decoupling takes place, i.e. 
mixing between the active and sterile neutrinos rapidly falls. 

We illustrated the enhancement of B(I/ in two scenarios: a model with one 

effective sterile neutrino that mimics the effect of a generic extensions of the SM includ¬ 
ing heavy sterile fermions, and in a minimal realization of the inverse seesaw scenario 
compatible with current observations. Our results for upper bounds on B(T ^ 

\y e {<p, J/fi), V’(2*S'), T(15'), T(25'), T(35'), Z}] are still considerably smaller than the cur- 
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Table 2: Upper bounds B(U —> Comparison of the results reported in the literature with the 

bounds obtained in this work by using two different models (the last two columns). The bounds for other 
similar decay modes that have not been discussed in the literature can be found in Table 

rent experimental bounds (when available), but that situation might change in the future 
as more experimental research will be conducted at Belle II, BESIII, LHC, and hopefully 
at FCC-ee (TLEP). If one of the decays studied here is observed and turns out to have a 
branching fraction larger than the upper bounds reported here, then sources of LEV other 
than those coming from mixing with heavy sterile neutrinos must be accounted for. 
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Appendix A: Wilson Coefficients 

In this Appendix we present detailed expressions for the Wilson coefficients. All computa¬ 
tions have been made in the Feynman gauge. Contributions coming from the penguin and 
self-energy diagrams are shown in Fig. whereas the box diagrams are shown in Fig. 


We use the standard notation, Xi = mj/rn^, Xt = rn^/m^, Xq = (f /m^ = my/m^, 
and write 




(25) 


where r e { 7 ,^, box}. The coefficients Cy^l related to 7 and the box contributions are 
diagonal, Cy^l = 5ijCy\, while those related to the Z penguins can also involve a cou¬ 
pling to two different neutrinos, since the 3x3 mixing matrix is no longer unitary. We 
therefore separate the diagonal and nondiagonal parts of the corresponding coefficient 
^VL ~ where the second term depends on the parameter Cij defined by 


C«= E 


a=e,/j,,r 


(26) 


which, in the presence of sterile neutrinos, is generally different from Sij. Furthermore, from 
the plots presented in the body of the present paper we see that the region of 1714^5 ^ my 
is particularly interesting because there occurs the enhancement of the LFV decay rate. 
For the sake of clarity we thus expand our expressions in Xg and present here only the 
dominant terms. We also neglected, in the denominators of the loop integrals, the external 
momenta since they are negligible with respect to heavy neutrino masses. Therefore, up 
to terms O^x"^), our results read: 

y y _ 1 , -43a:f + lOSxf 6(5a:j - 6 )x‘f loga:* - Six, + 16 

288 d{Xi - 1)4 ’ ^ 



-1-1- 12 xi — llx^ + lOxf log a:* 
647r^(a:j — 1)^ 


+ cos^ OwC^ldi), 


(28) 


c^tixi,x,) = a, 
di - 


^qCii 


{xj - 2) (3 {xj - 1) + 2{xi - A)xi loga:») 
1287r^(a:i — 1)^ 

- l){xi{xi{2xi — 47) -|- 25) -|- 14) -|- 6{xi{12xi 
11527r2(xi - 1)4 


13) -I- 2) log Xi 


(29) 
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Figure 6: Penguin and self-energy diagrams contributing the LFV decay in Feynman gauge. 
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Figure 7: Box diagrams contributing the LFV decay ^ in Feynman gauge. 
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( 31 ) 


Appendix B: Formulas and hadronic quantities 

In this Appendix we collect the expressions used to constrain the parameters of the models 
discussed in the present paper, as well as the values of the masses and decay constants 
used in our numerical analysis. In the expressions below we used the value of Gf = 
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= 1.166 X 10 ^ GeV as extracted from /r ^ In our scenarios, in which we 

extended the neutrino sector by adding heavy sterile neutrinos, the Fermi constant becomes 
Gp = \ Uei\'^\Unj\'^. For the models used in this paper, we checked to see that 

Gp = G/j, remains an excellent approximation. 


fi ey: We use the experimentally established upper bound B(/i ^ ey) < 5.7 x 
10 “^^, and the expression 


B(^ ^ ey) = 
G^{x) = 


V2G' 


F^W^W ^5 


ml\U;,U,4G,{x,)\\ 


1287r5F^ 

2x^ + 5x‘^ — X 

4(1 — xy 2(1 — xy 


3x^ 


log x , 


(32) 


to get one of the most significant constraints in this study. Notice that we use 
= 1 — and we kept the dominant contribution with X 4 . 


• W —>■ igi' ’- Combining the measured B(kF —>■ ep) = 0.1071(16) and B(kF —>■ up) = 
0.1063(15), with the expression 


B(1F ^ iaP) 


y/^GpUliY 


247rF 
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EM 
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) 2- 
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w 


.jj2 I 

(33) 


we further restrain the possible values of 771.4 while varying the mixing angles in the 
largest possible range. 

. Ar^,. = R^/RZ - 1: The ratio of the leptonic decay widths of a given meson P, 

Rp — F(P ePe)/T{P iip^) was recently shown to be quite restrictive on the 
possible values of 7714^5 and r] j^. The most significant constraints actually come 
from Ar^r = 0.004(4) and Arp = —0.004(3), and the corresponding formula reads. 


Arp = —1 + 


ml{ml - ml) 


2\2 

fji) i 


[mU^li+ml) - {ml. - m^)^] xR^{ml,ml.,ml) 


ml{ml m2)2 ^|;7^.|2 [ml{ml. + ttiJ) - {ml. - 777^)2] y-G{m],, ml., ml) 

(34) 


• Z pp : To saturate the experimental F(E ^ invisible) = 0.499(15) GeV, we sum 
over the kinematically available channels involving active and sterile neutrinos, 


r{z pp) = 
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GJ 
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Gf 


12 \/ 2 - 


Timz 


-XR^{ml,ml,ml)\Gij\ 
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o 2 2 2 r 
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(36) 
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[1 eee: We use the experimental upper bound B(/r —>■ eee) < 1 x 10 ^^0, and the 
expression 23 


expression 

p/ X GpiTi'^Tn 

= 6144^ t; 
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(36) 


60 


with the loop functions Fg“®, F^^, defined in 

Finally, the values of hadronic quantities not discussed in the body of the paper but 
used in our numerical analysis are listed in Table 0 


Quantity Value Ref. 

Quantity Value Ref. 

1.0195 GeV 

16 


U 241(18) MeV 

65 


mj/^ 3.0969 GeV 

16 


m 418(9) MeV 

28 


^ij}( 2 S) 3.6861 GeV 

16 


Ui2S)/m 0.713(16) 

16 


mr 9.460 GeV 

16 


fr 649(31) MeV 

30 


ttiT( 2 S') 10.023 GeV 

16 


fr( 2 S) 481(39) MeV 

30 


mT( 35 ) 10.355 GeV 

16 


/t( 35 ) 539(84) MeV 

31 



Table 3: Masses and decay constants used in numerical analysis. 


® Notice that the ratio of decay constants fip(2S)/fj/tp has been obtained from the corresponding 
(measured) electronic widths and the expression T{‘tpn e+e“) = {27171^^). 
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